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The quantum correction to the conductivity in disordered quantum wires with linear Rashba 
spin-orbit coupling is obtained. For quantum wires with spin-conserving boundary conditions, we 
find a crossover from weak antilocalization to weak localization as the wire width W is reduced 
using exact diagonalization of the Cooperon equation. This crossover is due to the dimensional 
dependence of the spin relaxation rate of conduction electrons, which becomes diminished, when 
the wire width W is smaller than the bulk spin precession length Lso- We thus confirm previous 
results for small wire width, W/L S o < 1 [S. Kettemann, Phys. Rev. Lett. 98, 176808(2007)], where 
only the transverse 0-modes of the Cooperon equation had been taken into account. We find that 
spin helix solutions become stable for arbitrary ratios of linear Rashba and Dresselhaus coupling in 
narrow wires. For wider wires, the spin relaxation rate is found to be not monotonous as function of 
wire width W: it becomes first enhanced for W on the order of the bulk spin precession length Lso 
before it becomes diminished for smaller wire widths. In addition, we find that the spin relaxation 
is smallest at the edge of the wire for wide wires. The effect of the Zeeman coupling to the magnetic 
field perpendicular to the 2D electron system (2DES) is studied and found to result in a modification 
of the magnetoconductivity: it shifts the crossover from weak antilocalization to weak localization to 
larger wire widths W c - When the transverse confinement potential of the quantum wire is smooth, 
the boundary conditions become rather adiabatic. Then, the spin relaxation rate is found to be 
enhanced as the wire width W is reduced. We find that only a spin polarized state retains a finite 
spin relaxation rate in such narrow wires. Thus, we conclude that the injection of polarized spins 
into nonmagnetic quantum wires should be favorable in wires with smooth confinement potential. 
Finally, in wires with tubular shape, corresponding to transverse periodic boundary conditions, we 
find no reduction of the spin relaxation rate. 

PACS numbers: 72.10.Fk, 72.15.Rn, 73.20.Fz 



I. INTRODUCTION 

Spintronic devices which rely on coherent spin preces- 
sion of conduction electrons 1,2 require a small spin relax- 
ation rate. As the electron momentum is randomized due 
to disorder, spin-orbit (SO) interaction is expected to re- 
sult not only in a spin precession but in randomization 
of the electron spin, the D'yakonov-Perel' spin relaxation 
with rate l/r s . 3 This spin relaxation is expected to vanish 
in narrow wires whose width W is of the order of Fermi 
wavelength Af, 4 ' 5 since the back scattering from impu- 
rities can in one-dimensional wires only reverse the SO 
field and thereby the spin precession. In this paper, we 
show, however, that 1/t s is already strongly reduced in 
much wider wires: as soon as the wire width W is smaller 
than bulk spin precession length L so , which is the length 
on which the electron spin precesses a full cycle. This ex- 
plains the reduction of the spin relaxation rate in quan- 
tum wires for widths exceeding both the elastic mean- 
free path l e and A^, as recently observed with optical 6 
as well as with weak localization measurements. 7-11 Since 
L so can be several [im and is not changed significantly as 
the wire width W is reduced, the reduction of spin relax- 



ation can be very useful for applications: the spin of con- 
duction electrons precesses coherently as it moves along 
the wire on length scale L so . It becomes randomized 
and relaxes on the longer length scale L S (W) = \fD e r s 
only \D e = vpl e j2 (vp, Fermi velocity) is the 2D dif- 
fusion constant]. Quantum interference of electrons in 
low-dimensional, disordered conductors is known to re- 
sult in corrections to the electrical conductivity Act. This 
quantum correction, the weak localization effect, is a very 
sensitive tool to study dephasing and symmetry-breaking 
mechanisms in conductors. 12-14 The entanglement of spin 
and charge by SO interaction reverses the effect of weak 
localization and thereby enhances the conductivity. This 
weak antilocalization effect was predicted by Hikami et 
at 15 for conductors with impurities of heavy elements. 
As conduction electrons scatter from such impurities, the 
SO interaction randomizes their spin. The resulting spin 
relaxation suppresses interference in spin triplet config- 
urations. Since the time-reversal operation changes not 
only the sign of momentum but also the sign of the spin, 
the interference in singlet configuration remains unaf- 
fected. Since singlet interference reduces the electron's 
return probability, it enhances the conductivity, which is 
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named the weak antilocalization effect. In weak magnetic 
fields, the singlet contributions are suppressed. Thereby, 
the conductivity is reduced and the magnetoconductivity 
becomes negative. The magnetoconductivity of wires is 
thus related to the magnitude of the spin relaxation rate. 
In Sec. II, we first derive the quantum corrections to the 
conductivity for wires with general bulk SO interaction 
and relate it to the Cooperon propagator. In Sec. Ill, we 
diagonalize the Cooperon for two-dimensional (2D) elec- 
tron systems with Rashba SO interaction. We compare 
the spectrum of the triplet Cooperon with the one of the 
spin-diffusion equation. In Sec. IV, we present the solu- 
tion of the Cooperon equation for a wire geometry. We 
review the solutions of the spin-diffusion equation in the 
wire geometry and compare the resulting spin relaxation 
rate with the one extracted from the Cooperon equation. 
Then we proceed to calculate the quantum corrections to 
the conductivity using the exact diagonalization of the 
Cooperon propagator. In the last part of this section, 
we consider two other kinds of boundary conditions. We 
calculate the spin relaxation rate in narrow wires with 
adiabatic boundaries, which arise in wires with smooth 
lateral confinement and regard also tubular wires. In Sec. 
V, we study the influence of the Zeeman coupling to a 
magnetic field perpendicular to the quantum well in a 
system with sharp boundaries and analyze how the mag- 
netoconductivity is modified. In Sec. VI, we draw the 
conclusions and compare with experimental results. In 
Appendix A, we give the derivation of the non-Abelian 
Neumann boundary conditions for the Cooperon propa- 
gator. In Appendix B, we show the connection between 
the effective vector potential As due to SO coupling and 
the spin relaxation tensor. In Appendix C, we give the 
exact quantum correction to the electrical conductivity 
in 2D. In Appendix D, we detail the diagonalization of 
the Cooperon propagator. In the following, we set h = 1. 



II. QUANTUM TRANSPORT CORRECTIONS 

If the host lattice of the electrons provides SO interac- 
tion, quantum corrections to the conductivity have to be 
calculated in the basis of eigenstates of the Hamiltonian 
with SO interaction 



H 



-L(p + eAf + F(x) - 1 7 <t (B + Bso(p)) , (1) 
2m e 2 



where m e is the effective electron mass. A is the vector 
potential due to the external magnetic field B. B^ Q — 
(B SOx , B SOy ) is the momentum dependent SO field, er 
is a vector, with components <T i; i = x,y,z, the Pauli 
matrices, 7 is the gyromagnetic ratio with 7 = gp, B with 
the effective g factor of the material, and /i B = e/2m e is 
the Bohr magneton constant. For example, the breaking 
of inversion symmetry in III-V semiconductors causes a 
SO interaction, which for quantum wells grown in the 



[001] direction is given by 16 

- i 7 B so ,D = otx(-e x p x + e y py) +Jd((1xPxpI ~ e y P y pl)- 

(2) 

Here, ot\ — Jd(p1) is the linear Dresselhaus param- 
eter, which measures the strength of the term linear 
in momenta p x ,Py in the plane of the 2DES. When 
(pi) ~ l/a 2 > hp (a is the thickness of the 2DES and 
kp is the Fermi wavenumber), that term exceeds the cu- 
bic Dresselhaus terms which have coupling strength jp). 
Asymmetric confinement of the 2DES yields the Rashba 
term which does not depend on the growth direction 



^7 b so.r = a 2 {e x p y - e y p x ), 



(3) 



with a-2 the Rashba parameter. 17,18 We consider the stan- 
dard white- noise model for the impurity potential, ^(x), 
which vanishes on average (V(x)) = 0, is uncorrelated, 
(y(x)F(x')) = <5(x - x')/27ti/t, and weak, e F r > 1. 
Here, v — m e /(2ii) is the average density of states per 
spin channel and r is the elastic scattering time. Go- 
ing to momentum (Q) and frequency (u>) representation, 
and summing up ladder diagrams, to take into account 
the diffusive motion of the conduction electrons, yield the 
quantum correction to the static conductivity as 15 



D, 
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(Q), (4) 



where a, /3 = ± are the spin indices, and the Cooperon 
propagator C is for epr 3> 1 (ep, Fermi energy), given 
by 



Cu:=E-E>(Q = P + P') = T 
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£',p'-q 
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where the impurity averaged electron propagator is given 
in the first Born approximation by 



<?I(P) 



1 



E-H (p) 



(6) 



and §p{p) is its complex conjugate, respectively. H$ 
is the Hamiltonian, Eq. (1), without disorder potential 
V. The impurity vertex (the cross) is given by 1/2tti/t. 
Impurity averaging products of Green's functions of the 
type (g R g R ) and (g A g A ) yield small corrections of order 
1/cft. Thus, the problem reduces to the calculation of 
the correlation function 



E,p + 
*~ 



E 



£',p'-q 



2lXVT ^-^ 

q 



g§Jp + <l)g A , a ,(p'- q ), (7) 



which simplifies for weak disorder ej?r > 1 to 
dfl 1 



2tt 1-irf 



(8) 
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where 



^ ^p'+q,^"' ^p— q.cr- (9) 



For diffusive wires, for which the elastic mean- free path 
l e is smaller than the wire width W, the integral is over 
all angles of velocity v on the Fermi surface. Using 

(p + eA) 2 1 m „ , ., 
p - q + eA 



S = \{<T + 0>), 

Q = p + p', 

we obtain to lowest order in Q, 

£ = -v(Q + 2eA + 2m e dS) + (Q + 2eA)aer' 



■i 7 ( < 7 , -«r)B. 



(10) 



Here, the SO couplings are combined in the matrix 

(ID 



-ai+^ D kl -Q' 2 



a 2 ot\ - "fDK x 

Thus, the Cooperon becomes 



i 

1 + zr(v(Q + 2eA + 2m e aS) + H a , + H 2 



,(12) 



where H a i = — (Q + 2eA)aer' and the Zeeman coupling 
to the external magnetic field yields 



1 



7(0-' - <x)B. 



(13) 



It follows that for weak disorder and without Zee- 
man coupling, the Cooperon depends only on the to- 
tal momentum Q and the total spin S. Expanding the 
Cooperon to second order in (Q + 2eA + 2m e aS) and 
performing the angular integral which is for 2D diffusion 
(elastic mean-free path l e smaller than wire width W) 
continuous from to 2tt and yields 



C(Q) 



1 



D e (Q + 2eA + 2eA s ) 2 + H. 



(14) 
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The effective vector potential due to SO interaction, 
As = m e aS/e ( where a = (a) denotes the matrix 
Eq. (11), as averaged over angle), couples to total spin 
vector S whose components are four by four matrices. 
The cubic Dresselhaus coupling is found to reduce the 
effect of the linear one to ot\ := a.\ — m e "/]jeF/2. Further- 
more, it gives rise to the spin relaxation term in Eq. (14), 



H JD =D e {mle FlD ) 2 (S 2 x + Sl). 



(15) 



In the representation of the singlet, \S = 0;m = 0) = 
(I U) - I lt>)/\/2 = | ^) and triplet states \S = l;m = 
0) = (I U) + I ;t»A/2 = | =t), \S = l;m = 1) = | TT 
),\S = l;m = —1) = | IX), C decouples into a singlet 
and a triplet sector. Thus, the quantum conductivity is 
a sum of singlet and triplet terms 
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D e (Q + 2eA) 5 



(S=l,m\ C(Q) | 5 = 1, 



(16) 



m=0,±l 



With the cutoffs due to dephasing 1/t v and elastic scat- 
tering 1 /t, we can integrate over all possible wave vectors 
Q in the 2D case analytically (Appendix C). 
In 2D, one can treat the magnetic field nonperturba- 
tively using the basis of Landau bands. 15 ' 19-23 In wires 
with widths smaller than cyclotron length kpl 2 B (lg, 
the magnetic length, defined by Big = 1/e), the Lan- 
dau basis is not suitable. There is another way to 
treat magnetic fields: quantum corrections are due to 
the interference between closed time-reversed paths. In 
magnetic fields, the electrons acquire a magnetic phase, 
which breaks time-reversal invariance. Averaging over 
all closed paths, one obtains a rate with which the mag- 
netic field breaks the time-reversal invariance, 1/tb- Like 
the dephasing rate 1/t v , it cuts off the divergence aris- 
ing from quantum corrections with small wave vectors 
Q 2 < l/D e TB- In 2D systems, tb is the time an elec- 
tron diffuses along a closed path enclosing one magnetic 
flux quantum, D e Ts = l%- I 11 wires of finite width W 
the area which the electron path encloses in a time tb 
is Wy/D e Ts ■ Requiring that this encloses one flux quan- 
tum gives 1/tb = D e e 2 W 2 B 2 /3. For arbitrary magnetic 
field, the relation 



1/tb = D e (2e) 2 B 2 (y 2 



(17) 
of 



with the expectation value of the square 
the transverse position (y 2 ), yields 1/tb = 
(l-l/(l + W 2 /3l 2 B ))D e /l 2 B . Thus, it is sufficient 
to diagonalize the Cooperon propagator as given by 
Eq. (14) without magnetic field, as we will do in the next 
chapters, and to add the magnetic rate 1/tb together 
with dephasing rate 1/t v to the denominator of C(Q) 
when calculating the conductivity correction, Eq. (16). 



III. THE COOPERON AND SPIN DIFFUSION 
IN 2D 

The Cooperon can be diagonalized analytically in 2D 
for pure Rashba coupling, a\ = 0, 7u = 0. For this case, 
we define the Cooperon Hamilton operator as 



C~ 



D, 



Q 2 + 2Q so (Q y S x — Q x S y ) + Q 2 so (S 2 + S 2 X ), 

(18) 
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Figure 1: (Color online) 2D spectrum of H c , Ki = Qi/Qso- 



and we define accordingly the spin-diffusion Hamiltonian 

^SD 



= d t s + D e H SD s, 



(25) 



where the matrix elements of the spin relaxation terms 
are given by 25 ' 26 (Appendix B) 



i 



I SI J 



= r 7 2 ((B so (k) 2 )% - (Sso(k) l Sso(k) J » . (26) 



For pure Rashba SO interaction, the spin-diffusion oper- 
ator H srt is in momentum representation 2 



with Q so = 2m e a2 = 2tt / L SCn where L so is the spin 
precession length. In the representation of the singlet 
| <=^) and triplet modes, {| TT), | = i) : | 11)} it becomes 



( Q 2 
o 
o 

V o 







) 2 

lso 



Q 2 

V^QsoQ- 
o 



o 

v / 2QsoQ+ 
2QL + Q 2 
V^QsoQ- 



V2Q 

^sc 



soQ+ 

Q 2 



(19) 

with Q± = Q y ± iQ x . Diagonalization yields the gapless 
singlet eigenvalues and the three triplet Cooperon eigen- 
values with a gap due to the SO coupling (see Fig. 1), 



E S {Q)/D e = 
E To (Q)/D e = 

E T± {Q)/D e = 




(20) 
(21) 



O 2 

16^,(22) 

^so 



where Eg denotes the singlet eigenvalue and Et , Et± 
the three triplet eigenvalues. Notice that the two min- 
ima of the lowest triplet eigenmode are shifted to Q = 
±(v / T5/4)Q so with a minimal eigenvalue of E/D e = 
(7/16)Q 2 . As we show in the following, this gap in 
the triplet modes is directly related to the D'yakonov- 
Perel' spin relaxation rate l/r s . We can get a better 
understanding of the spin relaxation induced by the SO 
coupling and impurity scattering by considering directly 
the spin-diffusion equation for the expectation value of 
the electron-spin vector 24 



B (r,t) = -<^t(r,t)o^(r,f)>, 



(23) 



where ip^ = (^t_ , tp_ ) is the two-component vector of the 
up (+), and down (-) spin fermionic creation operators 
and tp the two-component vector of annihilation opera- 
tors, respectively. In the presence of SO coupling, the 
spin-diffusion equation becomes for vf | V r s |<C 1/r, 

1 



= <9 f s + — s - D e V z s 

T. 



7 (B-2r((Vv F )B so (p))) x s 



(24) 







i^Qso^x 

-pi^ + k 2 i2Q so ky I , (27) 



i2Q 5Q k x i2Q so k 



y d c t s 



with 1/D e r s = Q 2 . In the 2D case, diagonalization 
yields the eigenvalues 



E (k) 
E±(k) 



1 



D p t s 



(28) 



k 2 + ^^±^^Jl + 16-^. (29) 



2 D e r s 2D e r s 



Qso 



Thus, we find that the spectrum of the spin-diffusion op- 
erator and the one of the triplet Cooperon Hamiltonian 
are identical in 2D (Ref. 28) as long as time-reversal sym- 
metry is not broken. This confirms that antilocalization 
in the presence of SO interaction, which has its cause in 
the suppression of the triplet modes in Eq. (16), is indeed 
a direct measure of the spin relaxation. Mathematically, 
there exists a unitary transformation 



H c — U CD H sr ,U cri , 



-— — 

1 

— -4= o 

V2 72 



(30) 



(31) 



with the according transformation between spin-density 
components Si and the triplet components of the 
Cooperon density s, 



1 



: S TT> 



-j={s x + is y ) = sn.. 



(32) 
(33) 
(34) 



This is a consequence of the fact that the four-component 
vector of charge density p = (p + +p_ )/2 and spin-density 
vector S are related to the density vector p with the 
four components (V'aV's) /V^, where a, (3 = ±, by a uni- 
tary transformation. The classical evolution of the four- 
component density vector p is by definition governed by 
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the diffusion operator, the diffuson. The diffuson is re- 
lated to the Cooperon in momentum space by substitut- 
ing Q —> p p' and the sum of the spins of the retarded 
and advanced parts, a and <x', by their difference. Using 
this substitution, Eq. (18) leads thus to the inverse of the 
diffuson propagator 



H d := 



D~ 



D, 



= Q 2 + 2Q so (Q y S x - Q x S y ) + Q 2 so {S 2 + S 2 



(35) 

with S = (er' — er)/2, which has the same spectrum as 
the Cooperon, as long as the time-reversal symmetry is 
not broken. In the representation of singlet and triplet 
modes the diffusion Hamiltonian becomes 



( 2QL + Q 2 


\-V2QsoQ- 



V^QsoQ- 
QL + Q 2 









Q 2 





-V^QsoQ^ 
o 
o 

Qlo + Q 2 



Comparing Eqs. (19) and Eq. (36), we see that diagonal- 
ization leads to Eqs. (20)-(22). 

It can be seen from Eqs. (28) and (29) that in the case of 
a homogeneous Rashba field, the spin-density has a finite 
decay rate. However, if we go beyond the pure Rashba 
system and include a linear Dresselhaus coupling, the 
first term in Eq. (2), we can find spin states which do not 
relax and are thus persistent. The spin relaxation tensor, 
Eq. (26), acquires nondiagonal elements and changes to 



1 



(ki 



4rk 2 F 



\a 2 




-aia 2 
ha 2 



(37) 



a 2 



a , 



with a = \/ a\ -Fof. For Q = and ot\ 
we find indeed a vanishing eigenvalue with a spin-density 
vector parallel to the spin-orbit field, s = so(l,l,0) T . 
Moreover there are two additional modes which do not 
decay in time but are inhomogeneous in space: the per- 
sistent spin helices, 29-33 



1 



S | -1 | sin( — (x-y) 



2tt 




g333 \\-y/2S 



Figure 2: (Color online) Persistent spin helix solution of the 
spin-diffusion equation for equal magnitude of linear Rashba 
and linear Dresselhaus coupling, Eq. (38). 
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(36) in the Q = case, the spin relaxation rate becomes 



i = 2 P 2 X 

t s a 



at 



-r + D e {m 2 e t F ^ D y 



(39) 



The condition for persistence is thus rather 6t\ = a 2 . 
This has been confirmed in a recent measurement (Ref. 
33). The existence of such long- living modes has an ef- 
fect on the quantum corrections to the conductivity. In 
this case, a± = a 2 — O-o, there is only weak localiza- 
tion in 2D. 35 ' 36 In the next sections we will make use 
of the equivalence of the triplet sector of the Cooperon 
propagator and the spin-diffusion propagator in quantum 
wires with appropriate boundary conditions and show 
how long-living modes may change the quantum correc- 
tions to the conductivity. 



IV. SOLUTION OF THE COOPERON 
EQUATION IN QUANTUM WIRES 



Quantum Wires with Spin-Conserving 
Boundaries 



1 



2tt 



S V2\ cos [—(x-y) , (38) 



^so 



(Fig. 2) and the linearly independent solution, obtained 
by interchanging cos and sin. Here, L S o = 7t/m e -\/2ao- 
One has to keep in mind that this solution is not an 
eigenstate anymore in a quantum wire. However, we will 
show that there exist also long persisting solutions in a 
quasi-lD case. 

It is worth to mention that in the case where cubic 
Dresselhaus coupling in Eq. (2) cannot be neglected, the 
strength of linear Dresselhaus coupling a\ is shifted 34 to 
Q'i = a.\ — TO e 7£>£i?/2, as mentioned in Sec. II, and, e.g., 



The conductivity of quantum wires with width W < 
Lip = yL\j\p is without SO interaction dominated by the 
transverse zero-mode Q y — 0. This yields the quasi-lD 
weak localization correction. 37 However, in the presence 
of SO interaction, setting simply Q y = is not correct. If 
we consider spin-conserving boundaries, rather one has to 
solve the Cooperon equation with the following modified 
boundary conditions as derived in Appendix A (Refs. 5, 
21): 



-— n • ( Vf [ 7 Bso(k) • S]) - td n ) C\ ds = 0, 



(40) 



(i 



where (...) denotes the average over the direction of vj? 
and k which we rewrite using Eq. (A9) for the given ge- 
ometry as 

(-id y + 2e(A s ) y )C (x, y = ±yj - 0, Wx, 

(41) 

where n is the unit vector normal to the boundary dS 
and x is the coordinate along the wire. The trans- 
verse zero-mode Q y = does not satisfy this condi- 
tion. Therefore, it is convenient to perform a non-Abelian 
gauge transformation, 21,24 so that the transformed prob- 
lem has Neumann boundary conditions, and the trans- 
formed Cooperon Hamiltonian can therefore be diago- 
nalized in zero-mode approximation for quantum wires. 
Since in quantum wires these boundary conditions apply 
only in the transverse direction, a transformation acting 

in the transverse direction is needed: C — > C — U A Cll\, 
with Ua = exp(i2e(Ag) !/ j/). Then, the boundary condi- 
tion simplifies to —id y C(x, y = ±W/2) = 0, Wx, and the 
Hamiltonian changes to 

H c = Q 2 - 2QsoQx[cos{Q so y)S y - sin(Q so y)S z ] 

+Q 2 so [cos 2 (Q so y)S 2 + sin 2 (Q so y)S 2 z 

- sin(Q so y) cos(Q so y)(S y S z + S z S y )} (42) 

= (Q + 2eA s ) 2 . (43) 

where the effective vector potential A s , as introduced in 
Eq.(14), 



A s = — aa = — 



-0-2 
a 2 




(44) 




is transformed to the effective vector potential A s after 
the transformation Ua has been applied to the Hamilto- 
nian 

A s = — a(y)S 

e 

_ 77^ / -a 2 cos(Q so y) -a 2 sm(Q so y) 
e 



which varies with the transverse coordinate y on the 
length scale of L so . Now, we can see already that for 
narrow wires W < L so , this vector potential varies lin- 
early with y, A s ~ —m e a2Qsoy/ e i like the vector po- 
tential of the external magnetic field B. Thus, it fol- 
lows, that for W < L so , the spin relaxation rate is 
l/ r « ~ Qloiv 2 ) ~ Qso^ 2 /^2, vanishing for small wire 
widths. As announced at the beginning, we thus see that 
the presence of boundaries diminishes the spin relaxation 
already at wire widths of the order of L so . If we include 
only pure Neumann boundaries to the Hamiltonian H C: 
i.e., using the wrong covariant derivative, this would not 
affect the absolute spin relaxation minimum and it would 
be equal to the nonzero one in the 2D case. We give a 
more precise answer in the following. 



B. Zero-Mode Approximation 

For W < L v , we can use the fact that the nth trans- 
verse nonzero-modes contribute terms to the conductiv- 
ity which are by a factor W/nL v smaller than the 0- 
mode term, with n a nonzero integer number. There- 
fore, it should be a good approximation to diagonalize 
the effective quasi-one-dimensional Cooperon propaga- 
tor, which is the transverse 0-mode expectation value of 
the transformed inverse Cooperon propagator, Eq. (43), 
Hid = (0 | H c \ 0). It is crucial to note that Hid con- 
tains additional terms, created by the non-Abelian trans- 
formation, which shows that taking just the transverse 
zero-mode approximation of the untransformed Eq. (18) 
would yield a different, incorrect result. We can now di- 
agonalize Hid and finally find the dispersion of quasi-lD 
triplet modes 




64-^(1 + c so (c so - 2)) 

^so , 



(46) 

where c so and t so are functions of the wire width W as 
given by 

2sm(Q so W72) sin(Q so WQ 
Cso = 1 7\ — —j j *so = 1 7\ — 777 — • ( 47 ) 



Inserting Eq. (46) into the expression for the quantum 
correction to the conductivity Eq. (16), taking into ac- 
count the magnetic field by inserting the magnetic rate 
1 /tb (W) and the finite temperature by inserting the de- 
phasing rate 1/t v (T), it remains to perform the integral 
over momentum Q Xl as has been done in Ref. 34. For 
Q so W < 1, the weak localization correction can then be 
written as 



Act 



Vh 



ir 



Vh> 



^/H v + B*(W)/4 ^/H v +B*(W)/4 + H s (W) 

(48) 



V Hw 



^H v + B*(W)/i + H s {W)/2 

in units of e 2 /2iT. We defined Hw = and the 

effective external magnetic field 

*W=(l-(l + S)~V (49) 



The spin relaxation field H S {W) is for QsoW < 1, 



Hs(W) = -(Q so W) 2 H s , 



(50) 
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suppressed in proportion to (W/L so ) 2 similar to B*(W), 
Eq. (49). Here, H s = l/4eD e r s , with l/r s = 2p 2 F a%r. As 
mentioned above, the analogy to the suppression of the 
effective magnetic field, Eq. (49), is expected, since the 
SO coupling enters the transformed Cooperon, Eq. (43), 
like an effective magnetic vector potential. 38 Cubic Dres- 
selhaus coupling, however, would give rise to an addi- 
tional spin relaxation term, Eqs. (15) and (B7), which 
has no analogy to a magnetic field and is therefore not 
suppressed in diffusive wires although it is width depen- 
dent due to presence of modified Neumann boundaries. 
When W is larger than SO length L so , coupling to higher 
transverse modes may become relevant even if W < L v 
is still satisfied, since the SO interaction may introduce 
coupling to higher transverse modes. 39 We will study 
these corrections by numerical exact diagonalization in 
the next section. One can expect that in ballistic wires, 
l e > W, the spin relaxation rate is suppressed in analogy 
to the flux cancellation effect, which yields the weaker 
rate, 1/t s (W) = {W/Cl e ){D e W 2 /12Lf), where C = 
10. 8. 40-42 Before we investigate the exact diagonaliza- 
tion in the pure Rashba case, we consider an anisotropic 
field with linear Rashba and Dresselhaus SO coupling to 
see which form the long persisting spin-diffusion modes 
have in narrow wires. Also, here, we can take advan- 
tage of the equivalence of Cooperon and spin-diffusion 
equation as far as time-reversal symmetry is not violated. 
We find three solutions whose spin relaxation rate decay 
proportional to W 2 for a.2 ^ &i and which are persis- 
tent for a2 = a.\- The first solution is s = so(ai,a2,0) T 
for Q x = which is aligned with the effective SO field 
B so (k) = — 2jk x (ai, a-i-, 0) T . In this case, we have 
according to Eq. (50) H 3 _ RD1 (W) = (l/12)(Q so W) 2 H s , 
with Ql Q = (2m e (a 2 — a 2 )) 2 /a 2 and l/r s = 2p 2 7 a 2 r, 
a = \/df + a| . As mentioned above by transforming 
the vector potential As, Eq. (45), this alignment occurs 
due to the constraint on the spin-dynamics imposed by 
the boundary condition as soon as the wire width W is 
smaller than the spin precession length L S o- In addition, 
we find two spin helix solutions in narrow wires, 

(51) 

and the linearly independent solution, obtained by in- 
terchanging cos and sin in Eq. (5f). The form of this 
long persisting spin helix depends therefore on the ra- 
tio of linear Rashba and linear Dresselhaus coupling 
strength, Fig. 3, and its spin relaxation rate is diminished 

a-S i? a ,RD2/3 = (l/2)-ff a ,RDl- 

C. Exact Diagonalization 

The exact diagonalization of the inverse Cooperon 
propagator, as obtained after the non-Abelian trans- 
formation, Eq. (43), is performed in the basis of trans- 




Figure 3: (Color online) Long persisting spin helix solution of 
the spin-diffusion equation in a quantum wire whose width W 
is smaller than the spin precession length L S o for varying ratio 
of linear Rashba ct2 = a sin if and linear Dresselhaus coupling, 
<ai = a cos </3, Eq. (51), for fixed a and L so = n/rn e a. 



verse standing waves, satisfying Neumann boundary con- 
ditions, {l/VW, VV/fFcos ((nn/W)(y - W/2))} with 

n £ N*, and the plane waves cxp(iQ x x) with momentum 
Q x along the wire. The results of this calculation for dif- 
ferent values of the dimensionless wire width Q so W are 
shown in Fig. 4. The numerical data points are attributed 
to the different branches of the eigenenergy dispersion by 
comparing their eigenvectors. For small Q so W, the re- 
sult is in accordance with the 0-mode approximation: For 
small wire widths W, the z-component of the total spin, 
S z , is a good quantum number, as can be seen by expand- 
ing Eq. (45) in Q so y. Thus, one can identify the low- 
est modes with the transverse zero-modes of the triplet 
modes corresponding to the eigenvalues of S z , m = 0, ±1, 
in the rotated spin axis frame, denoting them as E{ t0 >n= o} 
and Ei t ± n — \. The minimum of the E^ t0n=0 y mode is 
located at Q x = 0. The minimum of E{ t _ } is located 
at finite Q x > 0, in agreement with the 0-mode approx- 
imation. For larger Q so W, the modes mix with respect 
to the spin quantum number and the transverse quanti- 
zation modes. As a consequence, energy level crossings 
which are present at small wire widths are lifted at larger 
widths, since the mixing of spin and transverse quantiza- 
tion modes results in level repulsion being seen in Fig. 4 
as avoided level crossings. The branches, -E{to,o} an d 
E{t-,o}j evolve into two modes which become degenerate 
at large values of Q so W. These two modes are the only 
ones whose energy lies below the energy minimum which 
we obtained for the 2D modes, E/D e = (7/16)Q| , for 
a finite if^-interval around K x = 0. Therefore, we can 
identify these modes with edge modes which are created 
by the Neumann boundary conditions. We can confirm 
that these are edge modes by considering their spatial dis- 
tribution, shown in Fig. 5. Therefore, even in the limit 
of large widths W, we do get in addition to the spec- 
trum obtained for the 2D system with open boundary 
conditions case the edge modes, whose energy is low- 
ered as seen in Fig. 4. The presence of these edge states 
and the difference to the 2D system with open bound- 
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Figure 4: (Color online) Dispersion of the triplet Cooperon 
modes for different dimensionless wire units QsoW: (a) 
QsoW = 2, (b) QsoW = 8, (c) QsoW = 12, (d) Q so W = 30, 
plotted as function of K x = Q x /Qso- For QsoW » 3, £{to,o} 
and evolve into degenerate branches for large K x . 

(For QsoW — 30, not all high-energy branches are shown.) 



ary conditions can be seen in Appendix D in the non- 
diagonal elements which are proportional to the width 
times the functions Eqs. (D26) and (D27). Even in the 
limit of wide wires there are nondiagonal matrix elements 
which give a significant contribution which cannot be ne- 
glected. The modes above E/D e = (7/16)<3g are ex- 




h N 
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arb. units 

Figure 5: (Color online) Probability density of the Cooperon 
eigenmodes in the wire for QsoW/tt = 30. (a) 3D plot, (b) 
density plot for one of the two lowest branches, showing their 
edge mode character, (c) 3D plot and (d) density plot of the 
density of the third lowest mode, which shows bulk character. 



tended over the whole wire system and can thus be char- 
acterized as bulk-states, as seen in Figs. 5 (c) and (d). 43 
In Fig. 6, we compare the results which we obtained in 
the 0-mode approximation with the results of the ex- 
act diagonalization. We plot the absolute minima of the 
spectra as function of the dimensionless wire width pa- 
rameter QsoW/n = 2W/L S o- We confirm the parabolic 
suppression of the lowest eigenvalues for narrow wires 
~ W 2 /Ll Q , obtained earlier. 24,34 We note that the os- 
cillatory behavior of the triplet eigenvalues as function 
of W, obtained in the 0-mode approximation, 34 is dimin- 
ished according to the exact diagonalization. However, 
there remains a sharp maximum of E t Q at QsoW/tt « 1.2 
and a shallow maximum of £" t _ at QsoW/tt w 2.5. As 
noted above, the values of the energy minima of E t o and 
E t - at larger widths W are furthermore diminished as a 
result of the edge mode character of these modes. 
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Figure 6: (Color online) Absolute minima of the lowest eigen- 
modes -E{ t0j o}, ^{t-.o}! an d £ l { t+ } plotted as function of 
QsoW/n — 2W/L so - We note that the minimum of _E{ t _ j0 } 
is located at ±K X 7^ 0. For comparison, the solution of the 
zero-mode approximation Eto is shown. 



1. Comparison to Solution of Spin-Diffusion Equation in 
Quantum Wires 

As shown above, the spin-diffusion operator and the 
triplet Cooperon propagator have the same eigenvalue 
spectrum as soon as time-symmetry is not broken. 
Therefore, the minima of the spin-diffusion modes, which 
yield information on the spin relaxation rate, must be 
the same as the one of the triplet Cooperon propagator 
as plotted in Fig. 6. In Ref. 27, the value at K x = 0, 
with Ki = Qi/Q S o, has been plotted, as shown in Fig. 7. 
We note, however, that this does not correspond to the 




Figure 7: (Color online) Lowest eigenvalues at K x = plotted 
against Q S oW/n. For comparison, the global minimum of the 
Cooperon spectrum for QsoW > 9 is plotted, F3. Curves 
Fi[n] are given by 7/16 + ({n/(Q so W/TT))VT5/4) 2 , n G N. 
F2 shows the energy minimum of the 2D case, F2 = Fi[n = 
0]. Vertical dotted lines indicate the widths at which the 
lowest two branches degenerate at K x = 0. They are given 
by n/ (\/l5/4); consider that the wave vector for the minimum 
of the Et_ mode is (^15/4)Q S o- 



global minimum plotted in Fig. 6. The two lowest states 
exhibit two minima as can be seen in Fig. 4: one local 
at K x = and one global, which is for large Q so W at 
K x w 0.88. The first one is equal to the results given 




Figure 8: (Color online) The quantum conductivity cor- 
rection in units of 2e 2 /2tv as function of magnetic field B 
(scaled with bulk relaxation field H„), and the wire width 
W scaled with l/Qso for pure Rashba coupling and cutoffs 
1/ D e Qso T v = 0.08, l/D e Qs Q T = 4: Comparison of the zero- 
mode calculation (grid without shading) to the exact diag- 
onalization where the lowest 21 triplet branches and seven 
singlet branches were taken into account. 



by Ref. 27. For the WL correction to the conductivity, 
however, it is important to retain the global minimum, 
which is dominant in the integral over the longitudinal 
momenta. 



2. Magnetoconductivity 

Now, we can proceed to calculate the quantum correc- 
tions to the conductivity using the exact diagonalization 
of the Cooperon propagator. In Fig. 8, we show the re- 
sulting conductivity as function of magnetic field and as 
function of the wire width W. Here, we have included 
for all wire widths the lowest seven singlet modes and 
the lowest 21 triplet modes. We choose this number of 
modes so that we included sufficient modes to describe 
correctly the widest wires considered with QsoW = 10. 
Thus, for the considered low-energy cutoff, due to elec- 
tron dephasing rate 1/t v of \/D e Q^ T v = 0.08 and the 
high energy cutoff l/D e Ql Q T = 4 due to the elastic scat- 
tering rate, we estimate that seven singlet modes fall 
in this energy range. Since for every transverse mode 
there are one singlet and three triplet modes, we therefore 
have to include 21 triplet modes, accordingly. We note a 
change from positive to negative magnetoconductivity as 
the wire width becomes smaller than the spin precession 
length L so , in agreement with the results obtained within 
the 0-mode approximation, as reported earlier, 34 plotted 
for comparison in Fig. 8 (without shading). At the width, 
where the crossover occurs, there is a very weak magne- 
toconductance. This crossover width W c does depend on 
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the lower cutoff, provided by the temperature-dependent 
dephasing rate 1/t v . To estimate the dependence of W c 
on the dephasing rate, we have to analyze the contribu- 
tion of each term in the denominator of singlet and triplet 
terms of the Cooperon. A significant change should arise 
if 



-(W = W C ) = — (52) 



Assuming that this occurs for small wire widths, 
Q so W < 1, as confirmed for the parameters we used, 
we apply Eq. (50) to Eq. (52) and conclude that 



W c ~ — (53) 



If we calculate the crossover numerically in the 0-mode 
approximation we get the relation plotted in Fig. 10 
which coincides with Eq. (53). We note that the change 
from weak antilocalization to weak localization may oc- 
cur at a different width W c than the change of sign in 
the correction to the electrical conductivity Aa(B = 0) 
occurs, W wlj . However, we find that the ratio W c /W WIj 
is independent of the dephasing rate and the spin-orbit 
coupling strength Q so . 

Furthermore while there is quantitative agreement with 
the 0-mode approximation in the magnitude of the mag- 
netoconductivity for all magnetic fields for small wire 
widths W < L so , there is only qualitative agreement 
at larger wire widths. In particular, the total magnitude 
of the conductivity is reduced considerably in compar- 
ison with the 0-mode approximation. We can attribute 
this to the reduction of the energy of the lowest Cooperon 
triplet modes due to the emergence of edge modes, which 
is not taken into account when neglecting transversal spa- 
tial variations, as is done in the 0-mode approximation. 
Therefore, the 0-mode approximation overestimates the 
suppression of the triple modes, resulting in an overes- 
timate of the conductivity. Similarly, the magnetic field 
at which the magnetoconductivity changes its sign from 
negative to positive is already at a smaller magnetic field, 
as seen by the shift in the minimum of the conductivity 
towards smaller magnetic fields (Fig. 9), in comparison 
to the 0-mode approximation (unshaded) in Fig. 8. This 
is in accordance with experimental observations, which 
showed clear deviations from the 0-mode approximation 
for larger wire widths, with a stronger magnetic field de- 
pendence than obtained in 0-mode approximation. 7-10 
Note that the nonmonotonous behavior of the triplet 
modes as function of the wire width, seen in Fig. 4, cannot 
be resolved in the width dependence of the conductivity. 




Figure 9: (Color online) The relative magnetoconductivity 
Aa(B) — Aa(B — 0) in units of 2e 2 /2-7r, with the same pa- 
rameters and number of modes as in Fig. 8. 
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Figure 10: (Color online) Width of wire WQso at which there 
is a crossover from negative to positive magnetoconductivity 
as function of the lower cutoff 1/ ' D e Ql T v . 



D. Other Types of Boundary Conditions 

1. Adiabatic Boundary Conditions 

When the lateral confinement potential V is smooth 
compared to the SO splitting, that is, if \pd y V <C A so = 
2kpa2, where Xf is the Fermi wavelength, the bound- 
aries do not preserve the spin, s in 7^ s out , Eq. (40), since 
the spin may adiabatically evolve as the electron is scat- 
tered from such a smooth boundary. 44 If this applies, the 
potential is adiabatic and the spin of the scattered elec- 
tron stays parallel to the field B so as its momentum is 
changed. This leads to the boundary condition for the 
spin-density 27 

S x \y=±W/2 = 0, (54) 

Sy\y=±W/2 = 0, (55) 

dyS z \y = ± W /2 = 0. (56) 

We can transform this boundary condition to the one of 
the triplet Cooperon by using the unitary rotation be- 
tween the spin density in the Si representation and the 
triplet representation of the Cooperon, Sj, Eq. (32), which 
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leads to the boundary condition 



1 

71 



(— s x + is y )\y = ±w/2 = snly=±iv/2 



d 



yS z \y=±W/2 — dyS=}\y=±W/2 



~^={ s x + iSy)\y=±W/2 = Sll\y=±W/2 



0, (57) 
0, (58) 
0. (59) 



Now. if we require vanishing magnetization for the ID 
case, then the diagonalization is done straightforwardly 
as already calculated in Ref. 27. We use a ba- 
sis which satisfies the boundary conditions and there- 
fore consists of ~sin(gy)(l, 0, 0) T , ~cos(qy)(0, 1, 0) T , and 
~sin(<7y)(0, 0, 1) T , with q = nir/W, n G N*. However, 
looking at the spin-diffusion operator [Eq. (27)], we see 
immediately that if we set k to zero and use the fact 
that s x ^ y must vanish at the boundary and s z has to be 
constant for the chosen k, we receive a polarized mode. 
Although this mode is a trivial solution, it differs from all 
other due to the fact that it has a finite spin relaxation 
time as Q so W vanishes. For the choice of basis for diag- 
onalization, this means: We set qi = riiir/W for respec- 
tive §i therewith one state is described by {^1,^2,^3} = 
{?i,n+p 2 ,n+p 3 }, neW, p 2 G {-1, 0, 1, • • ■ }, p 3 eN. 
In the case {m, = {n,n,n} all branches di- 
verge with reference to the eigenvalues in the limes of 
QsoW —> 0, so that the spin relaxation time goes to 
zero for small wires. 27 In contrast, there is an additional 
branch in the case of P2 = — 1 , p 3 = which has a finite 
eigenvalue and therefore finite spin relaxation time for 
small wire widths 



E 



DeQsc 



K 



-2 



1 



32(Q so Wf 



0[(Q so W) 4 



(60) 

The smallest spin relaxation rate for vanishing Q so W, 
l/ r s,ir>, which is given for k x = 0, is found to be an 
eigenstate polarized in z direction which relaxes with the 
rate 1/t S; i_d = 2/r s . It shows compared with the other 
modes a monotonous behavior as function of Q so W . If 
we allow magnetization for the ID case, then the com- 
bination pi = — l,p 2 = leads to a valid solution. For 
wide wires the smallest absolute minimum is the 2D min- 
imum E/D e = (7/16)<5so> there are no edge modes. But 
already at a width of QsoW = n/\/E « 1.4 all modes 
except the z-polarized exceed the rate 1/t s _uj. 



2. Tubular Wires 

In tubular wires, such as carbon nanotubes, and InN 
nanowires in which only surface electrons conduct, and 
radial core-shell InO nanowires, 46 the tubular topology of 
the electron system can be taken into account by periodic 
boundary conditions. In the following, we focus on wires 
where the dominant SO coupling is of Rashba type. If 
one requires furthermore that this SO-coupling strength 
is uniform and the wire curvature can be neglected, 45 



the spectrum of the Cooperon propagator can be ob- 
tained by substituting in Eq. (22) the transverse momen- 
tum Q y by the quantized values Q y = nlix /W ', n is 
an integer, when W is the circumference of the tubular 
wire. Thus, the spin relaxation rate remains unchanged, 
1/t s = (7/16)D e <5s . If then a magnetic field perpen- 
dicular to the cylinder axis is applied as done in Ref. 
45, there remains a negative magnetoconductivity due to 
the weak antilocalization, which is enhanced due to the 
dimensional crossover from the 2D correction to the con- 
ductivity Eq. (C2) to the quasi-one-dimensional behavior 
of the quantum correction to the conductivity [Eq. (48)]. 
In tubular wires in which the circumference fulfills the 
quasi-one-dimensional condition W < L v , the weak lo- 
calization correction can then be written as 

Aa = 



w 



-2- 



B*(W)/4 y/Ef 



y/H v + B* (WO/4 + 7H s {W)/l& 



B*{W)/A + H S {W) 

(61) 



in units of e 2 /2n. As in Eq. (48), we defined Hw = 
1/AeW 2 , but the effective external magnetic field differs 
due to the different geometry: Assuming that W < Ib, 
we have 34 



-=# e (2 e ) 2 BV 

T~B 



1 (BWX 

2 l^rJ 



and the effective external magnetic field yields 



B*{W) = (2e) 



/ BW 



V 2tt 

= (2e)(Sr tub ) 2 



(62) 
(63) 

(64) 
(65) 



with the tube radius r tub . The spin relaxation field H s is 
H s = l/AeD e T s , with l/r s = 2p 2 7 a 2 r, or in terms of the 
effective Zeeman field Bgo, 



57 B so (e F ) 2 
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(66) 



Thus the geometrical aspect, (y 2 ) t ubo/(2/ 2 ) P ianar ~ 6.6, 
might resolve the difference between measured and cal- 
culated SO coupling strength in Ref. 45 where a planar 
geometry has been assumed to fit the data. This assump- 
tion leads in a tubular geometry to an underestimation 
of H S (W). The flux cancellation effect is as long as we 
are in the diffusive regime, l e W, negligible. 



V. MAGNETOCONDUCTIVITY WITH 
ZEEMAN SPLITTING 

In the following, we want to study if the Zeeman term, 
Eq. (13), is modifying the magnetoconductivity. Accord- 
ingly, we assume that the magnetic field is perpendicular 
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to the 2DES. Taking into account the Zeeman term to 
first order in the external magnetic field B = (0,0, B) T , 
the Cooperon is according to Eq. (14) given by 



C(Q) = 



1 



D e (Q + 2eA + 2eA s ) 2 + ih(tr' - <x)B ' 



(67) 



This is valid for magnetic fields jB <C 1/r. Due to the 
term proportional to (a 1 — <x), the singlet sector of the 
Cooperon mixes with the triplet one. We can find the 



eigenstates of C 



with the eigenvalues 1/A.j. Thus, 



the sum over all spin up and down combinations a/3, /3a 
in Eq. (4) for the conductance correction simplifies in the 
singlet-triplet representation to 



^2 C <x/3j3a - ^2 (~ 

a/3 j 



.(^|i)(i|f±) + (TT TT) 

\i)(i\^) + (lim\ll})*i- (68) 



A. 2DEG 

The coupling of the singlet to the triplet sector lifts the 
energy level crossings at K = ±l/\/2 of the singlet E$ 
and the triplet branch Et_ as can be seen in Fig. 11 for a 
nonvanishing Zeeman coupling. The spectrum, which is 
not positive definite anymore for all wave vectors, K = 
Q/Qso, is given by 
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Figure 11: (Color online) (a) Real and (b) imaginary parts 
of the spectrum of the 2D Cooperon with Zeeman term of 
strength g/j, B B / D e Qi Q = 0.4. E B ,2D,i (black), E b ,2d,2 (red 
dashed), -E b ,2d,3 (green), _E b ,2d,4 (blue dashed). Dashed ver- 
tical lines are located at K x = ±l/\/2, the wave vector where 
the triplet mode Et_ and the singlet mode Es are crossing 
each other (without loss of generality K y = 0). 



Weak Field 



B - 4Ki 



1, 



fi 



h = 3 V3 x /l08A-4 + /s - 18K; 



B = gfi B B/D e Ql Q1 



(72) 
(73) 
(74) 



Thus, there are spin states with the same real part of 
the Cooperon energy, so that they decay equally in time, 
but the imaginary part is different, so that they precess 
with different frequencies around the magnetic field axis. 
A significant change of the Cooperon spectrum appears 
when g/j, B B/D e exceeds Ql Q , as can be seen in Fig. 12(b). 
All states with a low decay rate do precess now, due to 
a finite imaginary value of their eigenvalue. Associated 
with this change is also a change of the dispersion of 
the real part of -E B ,2D,3 which changes for K y = from a 
nearly quadratic dispersion in K Xl ao + aiK^. for B < 1 to 



one which changes more slowly as ao + aiK : 
a 3 Kl for B > 1 [see Fig. 12 (a)]. 



2/3 



4/3 



In the case of a weak Zeeman field, B <C 1, the sin- 
glet and triplet sectors are still approximately separated. 
A finite B <C 1 lifts however the energy of the sin- 
glet mode to E Bt 2D,2( K = °)/ D eQ 2 S o = B 2 /2 + 0(B 4 ), 
thus the singlet mode attains a finite gap, correspond- 
ing to a finite relaxation rate. The absolute minimum of 
two of the triplet modes is also lifted by Eb,2D,2(K = 
±715/4) /D e Ql Q = 7/16 + (3/4)B 2 + C(S 4 ), while their 
value is independent of B at K = 0. In contrast, the min- 
imum of the triplet mode Eb,2D,3i which approaches Et + 
in the limit of no magnetic field (see Fig. 1) is diminished 
to E b ,2dA K = 0)/D e Q 2 so = 2- B 2 /2 + 0(B 4 ). So, in 
summary, a weak Zeeman field renders all four Cooperon 
modes gapfull and that gap can be interpreted as a finite 
relaxation rate or dephasing rate as the Zeeman coupling 
mixes all the spin states, breaking time-reversal invari- 
ance. 



13 



Strong Field 

If we expand the spectrum in 1 / B <C 1, we find that all 
modes have the same gap proportional to the strength of 
the SO coupling, D e Ql Q , while two modes attain a finite 
imaginary part with opposite sign 

E B . 2DA /D e Q 2 so = 1 + K 2 X , (75) 
E B , 2D , 2 /D e Q 2 so = 1 + K 2 X +0(1/B), (76) 
E B . 2D - 3/i /D e Q 2 so = 1 + KItiB + 0{1/B). (77) 

Thus, a strong Zeeman field polarizes the spins and leads 
to their precession. The SO interaction, which is too 
weak to flip the spins, merely results in a relaxation of 
all modes, corresponding to a dephasing of the spin pre- 
cession. 



B. Quantum Wire with Spin-Conserving Boundary 
Conditions 

In the following, we want to study if a Zeeman 
field modifies the magnetoconductivity and can shift the 
crossover from positive to negative Magnetoconductiv- 
ity as function of wire width W. We have seen that 
for appropriate parameters the critical width W c is small 
compared with L so . Therefore we stay in the 0-mode ap- 
proximation to get a better overview of the physics. To 
do so, we first analyze the spectrum. 
The modes with low decay rates are situated at K x = 
and K x w ±1 for small widths and small enough Zeeman 
field, B < 1, as can be seen in Fig. 13. For K x = 0, we 
have 

E^,o/D e Ql Q = & + B* (l + ( 7 8) 
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Figure 12: (Color online) (a) Real and (b) imaginary parts 
of the spectrum of the 2D Cooperon with Zeeman term of 
the strength g^, B B / D £ Q% Q = ... 2 in steps of 0.25 (w.l.o.g 
K y = 0). The B independent mode Et is not shown. 



and for K x = ±1, 

As in the 2D case, we have a mode which is independent 
of the Zeeman field and the spectrum is equal to E t o 
with the eigenvector (0, 1,0, 1) T . Using this spectrum, 
we estimate the correction to the static conductivity in 
the case of a magnetic field which we include by means 
of a Zeeman term together with an effective magnetic 
field appearing in the cutoff 1/tb as described in Sec. II. 
The g = g/8m e D e factor is used as a material-dependent 
parameter. In Fig. 14, we see that for large enough 
g factor, the system changes from positive magnetocon- 
ductivity — in the case without Zeeman field and a small- 
enough wire width — to negative magnetoconductivity at 
a finite Zeeman field for the same wire. Hence, the ratio 
Wc/Wwl changes and one has to be careful not to con- 
fuse the crossover defined by a change of the sign of the 
quantum correction, WL— >WAL, and the crossover in the 
magnetoconductivity. To give an idea how the crossover 
W c depends on g and the strength of the Zeeman field 
we analyze two different systems as plotted in Fig. 15: 
The first one, plot (a), shows the drop of W c in a sys- 
tem as just described where we have one magnetic field 
which we include with an orbital and a Zeeman part. For 
small g we have Q S oW c (cj) = Q S oW c (g = 0) — const g 2 , 
where const is about 1 in the considered parameter 
space. In the second system [Fig. 15(b)], we assume that 
we can change the orbital and the Zeeman field sepa- 
rately. The critical width is plotted against the Zeeman 
field. To calculate W c , we fix the Zeeman field to a cer- 
tain value, horizontal axis in plot (b), while we vary 
the effective field and calculate if negative or positive 
magnetoconductivity is present. For different Zeeman 
fields Bz/H s we get different W c . We see that W c is 
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Figure 13: (Color online) (a) Real and (b) imaginary parts 
of the spectrum of the Cooperon with Zeeman term of the 
strength j(i B B / D e Q% = 0. . .0.8 in steps of 0.1 in a finite 
wire of the width QsoW = 0.5. The B independent mode E t x> 
is not shown. 



shifted to larger widths as the Zeeman field is increased, 
Q so W c (B z /H s ) = Q so W c {B z = 0) + const (B z /H s ) 2 , 
where const is about 1 in the considered parameter space, 
while Act(1/t b = 0) (not plotted) is lowered as long 
as we assume small Zeeman fields. If we notice that 
B z mixes singlet and triplet states it is understood that 
there is no gapless singlet mode anymore and therefore 
A<t(1/t b = 0) must decrease for low Zeeman fields. 
To estimate <?, we take typical values for a GaAs/AlGaAs 
system and assume the electron density to be n s = 
1.11 x 10 11 cm -2 , the effective mass m e /m eo = 0.063, 
the Lande factor g = 0.75 and an elastic mean- free path 
of l e = 10 nm in a wire with Q so W = 1, corresponding 
to W = 1.2 //m, if we assume a Rashba spin-orbit cou- 
pling strength of oti = 5 meVA. We thus get g rj 0.1 and 
find that the Zeeman coupling due to the perpendicular 
magnetic field can have a measurable, albeit small effect 
on the magnetoconductance in GaAs/AlGaAs systems. 




B/Hs 

Figure 14: (Color online) The magnetoconductance Aa(B) in 
a magnetic field perpendicular to the quantum well, where its 
coupling via the Zeeman term is considered by exact diago- 
nalization while its effect on the orbital motion is considered 
effectively by the magnetic phase shifting rate 1/tb(W) for 
wire width QsoW = 1, dephasing rate 1/t> = Q.06D e Qs O , 
and elastic-scattering rate 1/r = 4D e Qso- The strength of 
the contribution of the Zeeman term is varied by the material- 
dependent factor g = gfi B H s /D e Qso in the range g — ... 1.5 
in steps of 0.5: The system changes from positive magneto- 
conductivity (g = 0, green) to negative one (g = 0.5 . . . 1.5, 
continuous decrease of the absolute minimum). 



VI. CONCLUSIONS 

In conclusion, in wires with spin-conserving bound- 
aries and a width W smaller than bulk spin precession 
length L so , the spin relaxation due to linear Rashba 
SO coupling is suppressed according to the spin relax- 
ation rate (l/r s )(W) = (tt 2 /3) (W/L so f (l/r s ), where 
1/t s = 2p 2 7 a 2 >T. The enhancement of spin relaxation 
length L s = -J D e r s (W) can be understood as follows: 
The area an electron covers by diffusion in time r s is 
WL S . This spin relaxation occurs if that area is equal 
L 2 so , 38 which yields l/L 2 s „ \/ DeTs „ (W/ L so ) 2 / L 2 so , 
in agreement with Eq. (50). For larger wire widths, the 
exact diagonalization reveals a nonmonotonic behavior 
of the spin relaxation as function of the wire width of 
the long-living eigenstates. The spin relaxation rate is 
first enhanced before it is suppressed as the widths W 
is decreased. The longest living modes are found to 
exist at the boundary of wide wires. Since we identi- 
fied a direct transformation from the spin-diffusion equa- 
tion to the Cooperon equation, we could show that these 
edge modes affect the conductivity: the 0-mode approx- 
imation overestimates the conductivity for larger wire 
widths Q so W > 1 since it does not take into account 
these edge modes. They add a larger contribution to 
the negative triplet term of the quantum correction than 
the bulk modes do, since they relax more slowly. This 
also results in a shift in the minimum of the conduc- 
tivity towards smaller magnetic fields in comparison to 
the 0-mode approximation. The reduction of spin relax- 
ation has recently been observed in optical measurements 
of rt-doped InGaAs quantum wires 6 and in transport 
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Figure 15: (Color online) (a) Change of crossover width W c 
with g factor: The magnetic field is included as an effective 
field 1/tb and in the Zeeman term. The strength of the con- 
tribution of the Zeeman term is varied by the material de- 
pendent factor g = gfi B H s /D e Q^ Q . (b) Change of crossover 
width W c with Zeeman field: To calculate W c , we fix the Zee- 
man field to a certain value, horizontal axis, while we vary the 
effective field independently and calculate if negative or pos- 
itive magnetoconductivity is present. For different Zeeman 
fields Bz/H s , we find thereby a different width W c . Here, we 
set g/8m e D e — 1. In (a) and (b), the cutoff due to dephasing 
is varied: l/D e Ql T v = 0.04,0.06,0.08 (lowest first). 



measurements. 7 ~ 10 > 47 Recently in Ref. 11, the enhance- 
ment of spin lifetime due to dimensional confinement in 
gated InGaAs wires with gate controlled SO coupling was 
reported. Reference 6 reports saturation of spin relax- 
ation in narrow wires, W <C L so , attributed to cubic 
Dresselhaus coupling. 34 The contribution of the linear 
and cubic Dresselhaus SO interaction to the spin relax- 
ation turns out to depend strongly on growth direction 
and will be studied in more detail in Ref. 48. Including 
both the linear Rashba and Dresselhaus SO coupling we 
have shown that there exist two long persisting spin he- 
lix solutions in narrow wires even for arbitrary strength 
of both SO coupling effects. This is in contrast to the 
2D case, where the condition a± — oli-, respectively in 
the case where the cubic Dresselhaus term cannot be ne- 
glected, ai — w e 7ep-/2 = CV2, is required to find persis- 
tent spin helices 29 ' 30 as it was measured recently (Ref. 
33). Regarding the type of boundary, we found that 
the injection of polarized spins into nonmagnetic mate- 



rial is favorable for wires with a smooth confinement, 
XpdyV < A so = 2kpct2- With such adiabatic boundary 
conditions, states which are polarized in z direction relax 
with a finite rate for wires with widths Q S oW <C 1, while 
the spin relaxation rate of all other states diverges in 
that limit. In tubular wires with periodic boundary con- 
ditions, the spin relaxation is found to remain constant 
as the wire circumference is reduced. Finally, by includ- 
ing the Zeeman coupling to the perpendicular magnetic 
field, we have shown that for spin-conserving boundary 
condition the critical wire width, W c , where the crossover 
from negative to positive magnetoconductivity occurs, 
depends not only on the dephasing rate but also depends 
on the g factor of the material. 
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Appendix A: Spin-Conserving Boundary 

In the following we set 7 = 1. In order to generate a fi- 
nite system, we need to specify the boundary conditions. 
These can be different for the spin and charge current. 
Here we derive the spin-conserving boundary conditions. 
Let us first recall the classical derivation of the diffusion 
current density in the wire at the position r. The cur- 
rent density j at position r can be related by continuity 
to all currents in its vicinity which are directed towards 
that position. Thus, j(r,£) = (vp(r — Ax)) where an 
angular average over all possible directions of the veloc- 
ity v is taken. Expanding in Ax = l e ~v/v and noting 
that (vp(r)) = 0, one gets j(r,t) = (v(-Ax)Vp(r)) = 
— (i>i^e/2)Vp(r) = —D e Vp(v). For the classical spin- 
diffusion current of spin component s;, as defined by 
j Si (r, t) — vsi(r,t). there is the complication that the 
spin keeps precessing as it moves from r — Ax to r, and 
that the SO field changes its direction with the direction 
of the electron velocity v. Therefore, the 0th order term 
in the expansion in Ax does not vanish, rather, we get 

j Si (r, t) = (vs$(t, t)) - D e Vsi(r, t), (Al) 

where is the part of the spin-density which evolved 
from the spin-density at r — Ax moving with velocity v 
and momentum k. Noting that the spin precession on 
ballistic scales t < t is governed by the Bloch equation 

^- = § x B so (k) - -§, (A2) 
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we find by integration of Eq. (A2) that after the scat- 
tering time r, the spin-density components are given by 
= — r[B so (k) xs]j so that we can rewrite the first 
term in Eq. (Al) yielding the total spin-diffusion current 

as 

j Si = -t{v f [B so (k) x B ] 4 ) - D e \7 Sl . (A3) 

In this section, we consider specular scattering from the 
boundary with the condition that the spin is conserved, 
so that the spin current density normal to the boundary 
must vanish 



Appendix B: Relaxation Tensor 

To connect the effective vector potential As with the 
spin relaxation tensor, we notice that f can be rewritten 
in the following way: 

i- = r((B so (k) 2 )<5 ?J - (B BO (k) i S BO (k) i )) i= i..3j=i..3 

(Bl) 

using Uod, Eq. (30), 



n-Uhw/2 =0, (A4) 

where n is the vector normal to the boundary. Noting the 
relation between the spin-diffusion equation in the Sj rep- 
resentation and the triplet components of the Cooperon 
density §, ({| Tt>, I =0, U»), Eq. (30), 

U C n{eijkBso,j)i=1..3,k=l..3Uc D 

= -i((si\B so ■ S|sfc))i = i.. 3j fc=i..3, (A5) 

where the matrix U CD is given by Eq. (31), we can 
thereby transform the boundary condition for the spin- 
diffusion current, Eq. (A4), to the triplet components of 
the Cooperon density Sj , 

= n-j Si | u=±w -/ 2 . (A6) 

Requiring also that the charge density is vanishing nor- 
mal to the transverse boundaries, which transforms into 
the condition — i<9 n /5| Surfaco = for the singlet component 
of the Cooperon density p, we finally get the boundary 
conditions for the Cooperon without external magnetic 
field, Eq. (40), 

^__L n . ( VF [B so (k) ■ S]) - id n ^j q surface = 0. (A7) 

The last expression can be rewritten using the effective 
vector potential As, Eq. (14), 



= rUU{B so (k) x )S 2 x + (Sso(k),)5. y 2 

+ (B so (k) x B so (k) y )(S x S y + S y S x )}U CD (B2) 
= r[/t D ((B so (k).S) 2 )[/ CD (B3) 

= 4^o((vFB so (k).S) 2 )(7 CD . (B4) 

V F 

Because 

4((v F [B so (k).S]) 2 ) = %{(y F [B so (k).S])) 2 (B5) 

Vp vp 

is true for linear Rashba and linear Dresselhaus SO cou- 
pling but, in general, false if cubic-in-k terms are included 
in the SO field, we have to write 

r((B so (k).S) 2 ) = ^((v F B so (k).S)) 2 + ct (B6) 

Vp 

so that we conclude 

i- = UUD e (2eA s ) 2 + ct)U CD (B7) 

with the separated cubic part ct = D e m 2 ,e F '-f^(S x + S 2 ). 
This reflects nothing but the fact that the effective SO 
Zeeman term in Eq. (1) can only be rewritten as a vector 
potential As when the SO coupling is linear in momen- 
tum. 



(n • 2eA s - id n ) C| Surface = 0. (A8) 



Appendix C: Weak Localization Correction in 2D 



In the case of Rashba and linear and cubic Dresselhaus 

SO coupling in (001) systems, we get In contrast to the case where we have a wire with a 

finite width, we can calculate the weak localization cor- 

D e n . / /rt n \ ow rection to the conductivity analytically in the 2D case. 

— 2eA s = -(v F (B so (k) • Sj) . J J . J 2 

t rhe cutoffs due to dephasmg c\ = l/D e Q* T v and elastic 

(_( QL1 _ 7o("ie"F) 2 \ —ct2 \ scattering c 2 = \/D e Q\ r determine whether we have a 

4 m(m e v F ) 2 \ positive or negative correction. Integrating over all possi- 

a2 ai i / ble wave vectors K = k/Q so in the case without bound- 

(A9) aries yields 
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2p 2 i rV^z / i i 

Act = - — \ dK(2irK) ' 



1 1 



E T+ {QsoK)/Q 2 so + a E T _ (Q S0 K)/Ql + d 



2e2 '4infi + ^ + ±mii 



2jt \ 2 V ci J 2 V 1 + ci 



VT6+C2+1 



arctan | ) - arctan . 

4 Vto+ c i / V Vts+^i J 1, 2 + Ci 



111 



I I -J ) + arctan ' ^ ? e+C2 1 



4 VtB+ci / V V / TO+ C i / 1, / 1 + Cl 

111 



(CI) 



TE+°i 2 \| + ci + c 2 - 2y^+ c 2 

). (C2) 



As an example, we choose parameters which have 
been used in the case of boundaries, l/D e Ql Q T v = 
0.08, l/D e Q 2 so T = 4: Acr/(2e 2 /27r) = -0.29. The ex- 
act calculation of wide wires (Q S oW > 1) approaches 
this limit as can be seen in Fig. 8. The weak localiza- 
tion correction in 2D as function of these parameters is 
plotted in Fig. 16. 
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Figure 16: (Color online) Weak localization correction in 2D 



in units of (2e /2tt) 
and C2 = 1/DeQsoT- 



The parameters are ci = l/D e Q 
Thick line indicates Act = 0. 



2 

so 1 <p 



The spin subspace is thus represented by 4 x 4 matrices, 
which we order starting with the singlet 5 = and then 
S = l,m = 1, m = 0, and m = — 1. Thus, we get 



(Q x ,n | H c | Q x ,n) = Q 2 so 







( A 

c 



\ D 











B n iF n D n 
-iF n C n iF n 



(D2) 

The calculation of the matrix elements yields (we set 

P = QsoW/n) 



(D3) 
(D4) 

(D5) 

(D6) 

(D7) 
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Appendix D: Exact Diagonalization 

We write the inverse Cooperon propagator, the Hamil- 
tonian H c , in the representation of the longitudinal mo- 
mentum Q x , the quantized transverse momentum with 
quantum number n £ N, and in the representation of 
singlet and triplet states with quantum numbers S, m, 
where we note that H c is diagonal in Q x , 



(Q x ,n,S,m | H c | Q x ,n',S',m') 



(Dl) 



and for n > 0: 
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3 ^2 
7+A x 
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sin(P7r) 



4(n + P)(n-P) Ptt 



(D9) 
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s'm(Pn) Calculating the matrix elements for n = 0, n' > 0, we get 
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\/2(n'-2P)(n' + 2P)' 

-1 + (-1)"') cos(Ptt) 
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And for n > 0, n' > 0, we get 
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with the functions 
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